Magnetoconductivity in the presence of Bychkov-Rashba spin-orbit interaction. 
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A closed-form analytic formula for the magnetoconductivity in the diffusive regime is derived in 
the presence of Bychkov-Rashba spin-orbit interaction in two dimensions. It is shown that at low 
fields B ^ Bso, where Bso is the characteristic field associated with spin precession, D'yakonov- 
Perel' mechanism leads to spin relaxation, while for B ^ Bso spin relaxation is suppressed and 
the resulting spin precession contributes a Berry phase-like spin phase to the magnetoconductivity. 
The relative simplicity of the formula greatly facilitates data fitting, allowing for the strength of the 
spin-orbit coupling to be easily extracted. 



The conductivity in classically weak magnetic fields 
shows signatures of quantum interference that is consid- 
erably affected by the presence of spin-orbit (SO) interac- 
tions. Knowing the functional dependence of the magne- 
toconductivity, Aa{B), provides a sensitive tool for the 
extraction of the strength of the SO coupling. It is thus 
of great practical interest to obtain a simple, analytic 
formula for Aa{B). 

In a system with SO interactions the dominant spin re- 
laxation mechanism is the D'yakonov-Perel' (DP) mech- 
anism [ij. It describes how momentum relaxation by 
impurity scattering can lead to spin relaxation via the 
SO interaction. The spin relaxation rate is given by 
I/tso ^ ^so'^l^^ where A^o is the spin-splitting and 
r is the mean free time. For the special case when the 
SO interaction is of the Bychkov-Rashba kind, i.e., 0,0 
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the spin splitting at the Fermi surface is Ago = ^masoVp, 
where ago is the SO coupling. In this case, l/r^o = 
D{2maso/Kf' , where D — VpT/2 is the diffusion constant 
in two dimensions (2D) and vp '^s the Fermi velocity. 

Earlier works treated the effect of the SO interac- 
tion exclusively in terms of the DP relaxation 0, 0j 
i.e., the SO scattering rate was introduced as a cutoff 
1/tso in the triplet part of the interference processes 
leading to the well known anti-localization effect. How- 
ever, it became apparent on rewriting the Hamiltonian in 
Eq. ^ as Hso = ^ {p — PsoT )^ that the SO interaction 
which now appears as a spin-dependent vector potential 
T = i(z X (?) with charge Pso — 2maso could give rise 



to spin-dependent Aharonov-Bohm [a, [1 131 ^'^d Berry 
phase-like 9] effects. 

In 2D, the effects of various kinds of SO interactions on 
the magnetoconductivity have been studied extensively 
in the low field regime Ib ^ I 0, Ull . In this limit 
the magnetic length Ib — \Jh/2eB is many times larger 
than the mean free path / — vpT, and therefore the field 
range B <C Btr — H/AeDr is adequately described by the 
diffusion approximation It has been shown that an 
analytical expression for Aa{B) can be derived for the 



Bychkov-Rashba and the Dresselhaus kinds of SO inter- 
actions 0,^3 1 t^^t ^'^ analytical solution is possible 
when both terms are present together These solu- 
tions are expressed in the form of a series over effective 
Landau Level indices that are then summed numerically. 

In this paper, these summations are done analytically, 
thereby providing a closed-form formula for Aa{B). Al- 
though only the Bychkov-Rashba case is worked out in 
detail here, the solution is the same in the case when 
only the linear Dresselhaus term is present, i.e., when 
the cubic term is absent - which is the case in low den- 
sity systems. Various approximate formulas used in the 
literature are derived as limiting cases of the formula de- 
rived here and their physics highlighted. 

As is well known, quantum correction to the conduc- 
tivity arises from the interference of time reversed trajec- 
tories 0|. The amplitude of this interference, C(r, r'), is 
called the Cooperon. Using the formalism first developed 
m Ref. [13 (for a detailed review, see Ref. Ul]), it can be 
shown that in the diffusive regime the Cooperon in the 
presence of the Bychkov-Rashba interaction satisfies the 
equation: 7iC(r,r') — S^{r ~ r')/27ri/r^, where 
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The parameter v — 'm/2Trh'^ is the density of states 
per spin and the rate is introduced to account 

for dephasing. The spin-matrix S = z x S, where 
S = ^ [<T^ + (T^) is the total spin of the interfering waves, 
i.e., the retarded and advanced waves in the particle- 
particle channel (Cooper channel) [l5i |. 

In terms of the dimensionlcss quantities, Ti. = TL/Iuob) 
and C = {2TivT^'hiOBt)C , where the "cyclotron" fre- 
quency ojb — AeDB /h, the Cooperon equation reduces to 
iiC = 52 (r). In the circular gauge, A(r) = {B/2){z x r), 
the Hamiltonian Ti. expressed in terms of the raising and 
lowering operators reads 0, 0| : 



n = {aa^}-i^/2bZ{a^S+~aS-)+b,o{S^~S^J+b^ , (3) 

where, the notation {aat} = ^{aa^ +a^a) is used. The op- 
erators and a raise and lower the Landau level index n, 
and iS*^ raise and lower the Sz values, respectively. The 
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dimensionless variables bg 
where 



Bxn — 



Bso/B and = B^/B 

h 



AcDtso 'ieh 



and B,„ ~ 



4:eDT^ 



(4) 



The term proportional to V^so in Eq. Q mixes the 
Landau levels with the spin triplet states. The singlet 
S = sector as seen from Eg. ^ does not mix. As a re- 
sult F for = reduces to |lflIliE3 = {aa^l + 
with eigenvalues Es{n) = n + 1/2 + 6^. Further progress 
can be made by noting that the triplet sector conve- 
niently decomposes into 3x3 blocks [lO|, around each 
n spanned by the three vectors |n — Sz) ® \Sz) where 
Sz — —1,0,1 for n > 0. These blocks, labeled 7i„, take 
the form: 



+ bso 



i\/2bson 

O 







"iy/2bson 

~^^/2bso{n+l) e„+i + bso 



iy/2bso{n+l) , (5) 



bip. (The n = term is treated 



where, e„ = n + 1/2 
separately). 7i„ is easily diagonalized by taking linear 
combinations |n, to) = J2s ^ri-s^ Sz I" ~ '^•^) ^ I'^z)- The 
corresponding eigenvalues are labeled Et.m{n)^ where 
TO = -1,0, 1. 

The conductivity correction, 5(j{B), defined as |lfil |: 



5a{B) 
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where, a and are spin indices, can be expressed in 
terms of the singlet Es{n) and triplet Et^m{n) eigenvalues 
as E3: 
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To further simplify, the sum of the inverse eigenval- 
ues in the triplet sector can be written as = 
I]m^tTm(") Hi ['fi-n]ii/\'Hn\ , whcrc \Hn\ is the de- 
terminant and [Tirijii sltc the minors of the diagonal ele- 
ments, giving 



3el + AbsoCr, 



ibsob^ 



n=l, 



{b'so 



l)e„ -f 2bl^{bso + 2b^) 



(8) 

Eq. (jHJ when substituted into Eq. {Tj) gives a series so- 
lution for the magnetoconductivity (with the n — Q term 



properly included) that was first obtained in Ref. [lollll| . 

In the following, the sum over n is performed analytically 
to give a closed form expression for Aa{B). It is worth 
noting that one arrives at the same equations if instead 
of the Bychkov-Rashba interaction the linear Dresselhaus 
term was present (see Ref. llOlV In the latter case the 
parameter bso takes the form |lOj bso = (fi./4eZ3)2ri^r, 
where fii = ^k{{k1) — fc^/4). Hence, as a function of bso 
the final results derived here are equally valid in both 



cases. 



The crucial step to do the n sum is to expand St as: 
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(Note that the n sum has been extended to include the 
ri = term. It cancels the apparent divergence that 
appears in the last term in Eq. 0.) The advantage of 
expanding St in this way is that since tn is linear in n, 
the sum over n can be done using the formula: 
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where, ip{z) is the di-gamma function, 7 is the Euler con- 
stant, and as is a field dependent cutoff parameter |l6| . 
The variables Us and Vs are easily obtained as [T^: 



/ 27r 
Vs = 2 (5 cos ie~ —{I - s) 

+ AbsoVs + (5bL + 46so6y - 1) 
where, S and 6 are defined as: 



(11a) 
(lib) 




(12a) 
(12b) 



The magnetoconductivity is defined as Af7(i?) = 
Sa{B) — S(t{0). To find the zero-field value, 6a{0), the 
sum over n in Eq. is replaced by an integral over mo- 
mentum q, with the replacement uJon ^ Dq^. This gives 
the final result, 
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where, cto = /2nh. The constant C satisfies the re- 
quirement A(t(0) — 0. Eq. 1(1^ . combined with Eqs. 
(|ll|l and p2l) . is the main rcsuh of this paper. It pro- 
vides an analytic expression for Acr(i?) in the presence 
of the Bychkov-Rashba SO interaction. The formula, as 
mentioned earlier, is also valid for the case of the lin- 
ear Dresselhaus interaction. It is important to bear in 
mind that Eq. H13|) is only valid in the diffusive regime 
B <C Btr- In this limit Btr does not appear explicitly. 

The different limiting cases (i) B <C Bso and {ii) B ^ 
Bso are studied below: For B <C Bso, it can be shown 
from Eqs. and that: 5 w (i/V3)6so and 6 w 

cos~^(— i3V3), and therefore Us ~ 1 — i(4/-\/7)s and Vg ~ 
6so(3s2 + isV? — 2)/2. (It is assumed that Bso '3> B^, 
valid at low temperatures). Hence, in the low field limit, 
Eq. (|13|l reduces to the familiar form, first derived by 
Hikami, Larkin and Nagaoka (HLN) ^5j: 
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other hand, the SO interaction appears as a gauge field 
Bef f — B — 2sBso, akin to a Berry phase-like spin phase 
(see Ref. 0), and the sign of Aa{B) > 0, restoring 
positive magnetoconductivity. (Note that the s = sec- 
tor is not affected by the SO interaction in this limit.) 
This crossover from negative to positive magnetoconduc- 
tivity around B « Bso is an unambiguous signature of 
the presence of DP mechanism. 

To summarize, a simple closed-form analytic expres- 
sion has been derived for the magnetoconductivity in the 
presence of SO interaction with linear splitting of either 
the Bychkov-Rashba or the Dresselhaus kind. The rela- 
tive simplicity of the formula greatly facilitates data fit- 
ting from which the strength of the spin-orbit coupling 
can be extracted 0| . Note that although 1 /tso depends 
on the transport properties of the electrons, the scale Bso 
in Eq. depends only on pso = 2maso [13 ■ Hence, 
using Bso as a fitting parameter directly gives the value 

of Pso- 

AP greatly benefitted from discussions with A. M. 
Finkel'stein, M. Khodas, M. Manfra, R. de Picciotto, S. 
Schmult and S. H. Simon. 



In the high field hmit B > Bso, 5 
6 w 7r/6, Us ~ I and Vs ~ {l~ 2bsobip] 
Eq. 1)1 3|l in this limit reduces to: 



« (1 - 2fo,o6v)/V3, 
s. Up to a constant 



Aa(B) 



s=±l 



B - 2sBso 



- 21n 



B^p 

"b 



(15) 

As seen from Eqs. (|14|l and \ib\ the characteristic 
field scale Bso of the SO interaction plays fundamen- 
tally different roles depending on the strength of B. 
For B <C Bso the SO interaction as seen in Eq. H14() 
leads to spin relaxation via the DP relaxation mecha- 
nism, thus Bso in this case appears as a cutoff. At 
low temperatures, B^p <C Bso, the dominant term is 
Acr(B) « -%l}{l/2 + B^/B) + \ii{B^/B) resulting in neg- 
ative magnetoconductivity 5]. For B ^ Bso, on the 
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